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Time-dependent current of the electrons excited in the conduction band after uhrafast interband 
photogeneration is studied theoreticaUy. The transient photocurrent is calculated for the nonlinear 
regime of response to a stationary electric field. The response demonstrates transient absolute 
negative conductivity when the electrons are excited slightly below the optical phonon energy, while 
the periodic oscillations of the electric current appear after formation of the streaming distribution. 
The quenching of these peculiarities by the elastic scattering of electrons is also considered. 
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The steady-state nonlinear response of electrons with 
a strongly anisotropic (streaming) distribution appears 
due to the cyclic motion of electrons accelerating in 
the passive region (e < Tiljo) and rapidly emitting op- 
tical phonons when penetrating into the active region 
(e > TiLOo). This transport regime takes place under the 
conditions v <^ f^^ <C t'o, see Sec. 14.2 in Ref. 1 and Sec. 
35 in Ref. 2. Here and below, t^E — Po/\e\E is the time 
of ballistic flight of an electron across the passive region, 
E the strength of the applied electric field, e the electron 
charge, v the elastic scattering rate, Vo the rate of spon- 
taneous emission of optical phonons, po = y/2mhiOo the 
electron momentum corresponding to the optical phonon 
energy Tiujo, and m the electron effective mass. Stimu- 
lated THz emission of bulk holes in the streaming regime 
remains under investigation during the past 20 years, see 
Ref. 3 and references therein. However, a direct observa- 
tion of the current oscillations caused by the cyclic mo- 
tion of carriers is not possible for the steady-state case. In 
this paper we show that a transient streaming-oscillating 
(SO) photocurrent can be excited by an ultrashort inter- 
band pump creating electrons in the passive region after 
emission of the optical phonon cascade.'' 

Apart from the SO current, we have found that a tran- 
sient absolute negative conductivity, when the current 
flows in the direction opposite to the applied field, ap- 
pears if the excitation energy is close to the edge of the 
active region. The phenomenon of absolute negative con- 
ductivity (ANC) of photoexcited electrons in strong elec- 
tric fields has been discussed theoretically four decades 
ago^'^ for the case of stationary excitation. Now it is clear 
that the negative photoconductivity in this case cannot 
be observed neither in bulk materials^ nor in heterostruc- 
tures, because of accumulation of low-energy electrons 
with time and scattering of high-energy electrons by the 
low-energy ones,^ which leads to a fast broadening of 
the narrow energy distribution of photoexcited electrons. 
The case of ultrafast interband photoexcitation is free 
from these disadvantages because the accumulation takes 
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FIG. 1: Temporal evolution of electrons excited at t = with 
£ex ~ 0.64 htOo- (a) Position of ID electrons in the passive 
region at different instants in the case of (5-shaped initial dis- 
tribution, (b) Corresponding transient current showing the 
regimes of negative and oscillating conductivity. The dashed 
line is the current in the case of a Gaussian-shaped, with the 
width 0.2 ftuo, initial distribution, (c) Evolution of the initial 
(5-shaped distribution of 2D electrons: dotted circle at t = 0, 
dashed circle at t = ti, and solid segment plus a short straight 
line near the center (streaming component) at t > ti. 



place only at long times associated with energy relax- 
ation, and the electron-electron scattering remains inef- 
fective if the density of excited electrons is low. In con- 
trast to the steady-state negative conductivity, the effect 
we condider is transient, it exists in a finite time inter- 
val after the photoexcitation. Nevertheless, the origin of 
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this effect is the same as in the steady-state case, and the 
necessary condition for its reahzation is the partial inver- 
sion of the distribution of excited electrons. Recently, a 
transient ANC and negative absorption of photoexcited 
electrons have been suggested® for the regime of linear 
response, so the consideration of the nonlinear transient 
response is now timely. 

To describe the mechanisms of the transient ANC at 
t < tsE and oscillating current att > fes, we consider the 
model of one-dimensional (ID) conductor [Fig. 1 (a,b)], 
where electrons are excited at t = with the momenta 
ipex ~ ^V^TTT^ex Corresponding to the excitation energy 
Sex < Tt^o- Let us assume that p^x > Po/'^ and take into 
account that the electric field accelerates the electrons 
according to the dynamical equation pt = Pt=o + \e\Et 
(the electric field in Fig. 1 is assumed to be applied in 
the negative direction so that the electrons move from 
the left to the right). The electrons of the right group 
(with pt=o = Pex) come to the edge of the active re- 
gion at ti = {po — Pex)/\e\E and immediately drop to 
the point p = 0, where they do not contribute to the 
current. At the same instant, the electrons of the left 
group (with pt=o = ^Pex) have negative momentum, 
Pt=ti = Po — '^Pex ) so that the current abruptly becomes 
negative. Since the momentum of the electrons of the 
left group increases again, the current changes its sign 
from negative to positive at to = ^fi_E/2, when the abso- 
lute values of the momenta for the two groups become 
equal. The interval of the ANC is equal to to — ti and 
becomes maximal (equal to tgE/'^) when ti 0, i.e. 
when the electrons are cxcitcid very close to the edge. At 
t > Pex/\e\E both groups of electrons undergo cyclic mo- 
tion along the right-hand branch of the energy parabole, 
with the period Iqe- This leads to periodic oscillations of 
the current with the same period [Fig. 1 (b)]. The cur- 
rent, expressed in the figures here and below in units of 
jo = \e\nexPo/''n, where riex is the excited electron den- 
sity, linearly increases with time and drops abruptly each 
instant when one of the groups reaches the boundary of 
the active region (this occurs at t = (kpo T Pe2;)/|e|-E, 
where k is integer). If p^x < Po/'^, the oscillations re- 
main but the current is always positive. The main effects 
discussed above remain valid when the initial electron 
energy distribution is a broad peak, though the sharp 
features of the current become smoothed. 

In the 2D and 3D cases (which differ from each other 
only by geometrical factors), the evolution is more com- 
plicated because the processes of optical plionon emis- 
sion returning electrons to the point p = persist in a 
wide time interval, ti < t < ti + 2pex/\e\E for the 6- 
shaped initial distribution fe"^"^ cx 6{£ — Sex)- This inter- 
val begins when the surface jptj = Pex (here and below, 
Pt = p — eEt) touches the boundary of the active re- 
gion, IpI = Po, [sec Fig. 1 (c)] and ends when all the 
excited electrons emit optical phonons. Starting from ti , 
the distribution in the passive region acquires a continu- 
ous streaming part (at t > ti -I- 2pex/\e\E only this part 
remains). Nevertheless, both the ANC and the transient 



SO current exist in the 2D and 3D cases, and their prop- 
erties are similar to those discussed for the ID case. 

To describe the electron distribution in the passive re- 
gion, p < Po, we first analytically consider the case of 
ballistic transport over the passive region and then nu- 
merically calculate the current in the presence of elastic 
scattering described by the coUision integral Jeiif\pt)- 
The distribution function fpt is governed by the kinetic 
equation 
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Ut = Gt5{p) + JelU\pt) (1) 



with the initial condition /pt=o = fe'^^\ where e = 
£p = p^ /2ra,^^ and the normalization condition riex — 
2 dpfpt/{2TTnY, where d =2 or 3 for the 2D or 3D 
case and Vo is the area or volume of the passive region. 
Equation (1) assumes that the electrons instantcniously 
emit optical phonons so that their penetration into the 
active region is neglected and the generation term on the 
right-hand side, Gt, describes the electrons coming from 
the boundary of the active region. So, directly to the 
point p = 0. Using the conservation of the total number 
of electrons in the processes of optical phonon emission 
and elastic scattering, we find = eE • Jg dsfpt, where 
the integral over the boundary So appears owing to the 
relation /^^ dpVp/pt = J^^ dsfpt. 

Even in the ballistic regime, the presence of the inte- 
gral over So in the generation term makes Eq. (1) an 
integro- differential equation. One can solve this equa- 
tion analytically by separating the streaming contribu- 
tion to the electron distribution function according to 
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.gpBiit<5(Pfi±), where P6_L and p6|[ are the 
components of the momentum perpendicular and an- 
tiparallel to the electric field, respectively. The first part 
of this function satisfies the kinetic equation (1) without 
the right-hand side, while the second part is governed by 
the differential equation 
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where e is the unit vector in the direction of E and the 
initial condition is 3pg||t=o = 0. Using the trajectory 
method, one may transform Eq. (2) into an algebraic 
finite-difference equation which can be solved at p^ = Po- 
Thus, the streaming contribution takes the form 
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where km is the integer part of (\e\Et — p^) /po- De- 
pending on the variables P6|| and t, the sum in Eq. (3) 
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contains just one or two terms, because f^"^^ is nonzero 
only in the region e < hcOg. 

The current density, = 2(e/m) Jy rfpp/pt/(27r?i)'*, 
is written as jt = ejt with 
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where wc have used Eq. (3) and denoted Pm = 
min(po, |e|M). Analytical integration in Eq. (4) is 
possible for the ^-shaped initial distribution, /i^^' = 
{nexl PdD{eex)\^{e - Eex), whcrc pdoisex) is the d- 
dimensional density of states at the excitation energy. If 
\e\Et < Po — Pex, the current linearly increases with time, 
jt = e'^riexEt/m, while at \e\Et > Po — Pex its temporal 
dependence becomes complicated. 




FIG. 2: Transient current in 3D (a,b) and 2D (c,d) conductors 
for 5-shaped initial photoexcited distribution (solid) and for 
a Gaussian peak fJ'^'^' oc exp[— (e — £e2^)^/A^] (see Ref. 10) 
with the width A =0.2 TwJo (dashed). The excitation energies 
are Sex = 0.64 hu)o (a,c) and 0.9 hu)o (b,d). 

The results of numerical calculations of the current ac- 
cording to Eq. (4) for the 3D and 2D cases are shown 
in Fig. 2. We have chosen two excitation energies, close 
to the edge of the active region and far from this edge, 
and plot the temporal dependence of the current both for 
infinitely narrow (5-shaped) distribution of excited elec- 
trons and for a wide, Gaussian-shaped distribution. The 
transient current at large t oscillates around the value of 



the steady-state streaming current, jo/2. The plots for 
the 5-shaped excitation show sharp (kink-like) features 
at t = {kpo ^pex)/\£\E, whose origin is explained in the 
above discussion of the ID case (Fig. 1). The plots for 
the wide excitation close to the edge show similar sharp 
features dX t = kte,E, which are caused by the abrupt 
drop of the electron distribution at p = Po (these fea- 
tures are absent if the electrons are excited far from the 
edge). For the excitation close to the edge, there exist 
regions of ANC. One can find a simple criterion describ- 
ing the existence of ANC at small times: jt^o < at 

fhJa > nex[dhuJoPdD{lT-(^o)]^^ , where =2 or 3 is the 
dimensionality of the electron gas. In other words, the 
ANC at small times always takes place if the number of 
electrons excited close to the edge of the active region 
is sufficiently large. Comparing the plots (b) and (d) in 
Fig. 2, one can see that the absolute value of the negative 
response is larger for 2D electrons. 

The elastic scattering tends to isotropisc the distribu- 
tion function of electrons. For this reason, it is important 
to study the effect of this scattering on both the ANC and 
SO current. To do this, we consider Eq. (1) with the col- 
lision integral on the right-hand side. We restrict ourself 
by the 2D case, when the energy-independent scattering 
rate, v, can be used and this integral acquires the most 



simple form Jeii.f\pt) = -vif-pL - /if), where fH' (with 
e = Ep = p^ /2m) is the isotropic part of /pt, which is 
obtained by averaging /pt over the angle of the vector p. 
Then Eq. (I) is reduced to the integral equation 
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which can be transformed to a closed integral equation 

(s) 

for /e( • We have solved this equation numerically, by 
iterations. The results presented in Fig. 3 (a) show that 
the scattering suppresses both the magnitude and time 
interval of the ANC, while Fig. 3 (b) demonstrates how 
the amplitude of the oscillations decreases with time and 
with increasing scattering. An essential point is that the 
suppression of the ANC is not so strong as one may ex- 
pect, especially as concerns the time interval. If the scat- 
tering rate is v =50 t^^, the ANC interval decreases only 
by half in comparison to the ballistic case and becomes 
approximately 8 times greater than v^^. However, the 
scattering considerably suppresses the absolute value of 
the current. On the other hand, the suppression of the 
SO current is already strong at ~ f^^, and the oscilla- 
tions become washed out as p increases. 

Let us discuss the results in terms of numerical pa- 
rameters. The ballistic expression (4) is suitable for 
describing the transient response at the times smaller 
than the collision times. In clean samples, it is limited 
by the acoustic-phonon scattering rate Vac, according to 
t < ■ The measured low-temperature collision times 
for 2D electrons in clean GaAs-based structures'^ are of 
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FIG. 3: Transient current in 2D conductors in the presence 
of elastic scattering in the ANC (a) and SO (b) temporal 
regions. The initially excited distribution is the Gaussian of 
width A =0.1 TiuJo centered at e^^ =0.9 TvjJo- The dashed 
curves show the ballistic case (i^ = 0), while the solid curves 
are marked with the corresponding values of i^tes. 



about 150 ps. As follows from Fig. 2 (d), the region of 
ANC terminates at i ~ 0.35 t^E- Therefore, the influence 
of the scattering on the ANC is not essential at £^ > 0.35 
VacPo|\p\^ which givcs us an estimate i? > 4 V/cm for 
GaAs quantum wells. The numerical calculations accord- 
ing to Eq. (5) show that the ANC can be as well observed 
in the fields smaller by one order of magnitude, since the 
elastic scattering does not lead to a strong suppression of 
the effect; the time interval where the ANC exists in this 
case extends over a nanosecond scale. In order to gen- 
erate several periods of the transient oscillations in the 



streaming redime, the electric field should be consider- 
ably higher, i? > 10 V/cm. In bulk semiconductors,^ 
at low temperatures exceeds 100 ps (moreover, for InSb 
> 1 ns) so that the estimates given above remain 
valid. Since the frequency of the SO current, 2Tr/te,E, 
is proportional to E, one can obtain high-frequency os- 
cillations by increasing E up to Poi^o/\s\- In particular, 
application of high fields, E ~ 3x10^ V/cm, to wide-gap 
materials such as GaN (where po and z/q are large) can 
generate oscillations of about 0.5 THz, while for GaAs 
the attainable frequencies are several times smaller. 

Now we discuss the main assumptions. The neglect 
of inelastic scattering in the passive region assumes that 
the energy given to electrons via their acceleration by the 
field is much higher than the energy which they emit or 
absorb when scattering by acoustic phonons. This ap- 
proximation is well justified even at t ^ because 
of smallness of the sound velocity s and remains valid 
under the condition E ^ Vacfns/\e\. The neglect of a 
finite energy of electron penetration into the active re- 
gion is justified if this energy is small in comparison to 
hujo- This requires i^otsE ^ 1, which is easily attain- 
able since the typical times of optical phonon emission 
in quantum wells and bulk materials do not exceed 0.4 
ps. Next, since we consider the interaction of electrons 
with dispersionless optical phonons, the presented results 
have a limited applicability for heterostructures. Finally, 
the neglect of electron-electron scattering can be justi- 
fied if the excited electron density Uex is small enough so 
that this scattering mechanism is less important than the 
elastic scattering. Comparing the elastic scattering time 
of the order of 100 ps to the electron-electron scattering 
time,^ one can find Uex < 10^^ cm^'' for bulk GaAs and 
Hex < 10^ cm~^ for GaAs quantum wells. 

In conclusion, we have proposed and theoretically ver- 
ified a mechanism of the transient absolute negative con- 
ductivity, which is related to the process of transforma- 
tion of the nearly isotropic initial electron distribution 
to the strongly anisotropic (streaming) one. This effect 
appears to be robust against the elastic scattering, can 
exist over one nanosecond, and can be investigated by 
standard experimental methods. We have also described 
transient oscillations of the electric current caused by the 
streaming effect. In strong electric fields, the frequencies 
of these oscillations belong to the THz spectral region. 



[1] V.F. Gantmakher and Y.B. Levinson, Carrier Scattering 

in Metals and Semiconductors (Elsevier Sci. Ltd., 1987). 
[2] F.T. Vasko and O.E. Raichcv, Quantum Kinetic Theory 

and Applications (Springer, N.Y., 2005). 
[3] J. Shah, Ultrafast Spectroscopy of Semiconductors and 

Semiconductor Nanostructures, (Springer, N.Y., 1996). 
[4] M.A. Odnoblyudov, I.N. Yassievich, M.S. Kagan, and 

K.A. Chao, Phys. Rev. B 62, 15291 (2000). 
[5] H.J. Stocker, Phys. Rev. Lett. 18, 1197 (1967). 
[6] V.F. Elesin and E.A. Manykin, JETP Lett. 3, 15 (1966); 



Sov. Phys. - JETP 23, 917 (1966). 
[7] P.G. Harper, J.W. Hodby, and R.A. StradUng, Rep. Prog. 

Phys. 36, 1 (1973). 
[8] S.E. Esipov and LB. Levinson, Sov. Phys. - JETP 59, 

1113 (1984); Adv. Phys. 36, 331 (1987). 
[9] F.T. Vasko, JETP Lett. 79, 431 (2004); F.T. Vasko, A. 

Hernandez-Cabrera, and P. Aceituno, Phys. Rev. B 71, 

245325 (2005). 

[10] If the initial excitation creates electrons both in the pas- 
sive and active regions, with the distribution function 



5 



F'r\ wc have /i^-' = e{Tu^„ - e)Er=o^SL„' ^^crc 
the terms with fc 7^ describe the electrons rapidly trans- 
ferred to the passive region via optical phonon emission. 
Instead of using the initial condition, one may add a gen- 
eration term 5{t)fe''J^^ describing instantenious filling of 



the passive region to the right-hand side of Eq. (1). 
[11] P.J. Burke, LB. Spielman, J. P. Eisenstein, L.N. Pfeiffer, 
and K.W. West, Appl. Phys. Lett. 76, 745 (2000). 



